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Abstract 



O 

. We have experimentally demonstrated the interferometric complementarity, 



which relates the distinguishability D quantifying the amount of which-way 



(WW) information to the fringe visibility V characterizing the wave feature 

& . 

of a quantum entity, in a bulk ensemble by Nuclear Magnetic Resonance 
(NMR) techniques. We primarily concern on the intermediate cases: partial 
fringe visibility and incomplete WW information. We propose a quantitative 
measure of D by an alternative geometric strategy and investigate the relation 
between D and entanglement. By measuring D and V independently, it turns 
out that the duality relation D 2 + V 2 = 1 holds for pure quantum states of 
the markers. 
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I. INTRODUCTION 



Bohr complementarity [1] expresses the fact that quantum systems possess properties 
that are equally real but mutually exclusive. This is often illustrated by means of Young's 
two-slit interference experiment, where "the observation of an interference pattern and the 
acquisition of which- way (WW) information are mutually exclusive" [2]. As stated by Feyn- 
man, the two-slit experiment "has in it the heart of quantum mechanics. In reality it 
contains the only mystery" [3]. Complementarity is often superficially identified with the 
'wave-particle duality of matter'. As its tight association with the interference experiment, 
the terms of the "interferometric duality" or "interferometric complementarity" are more 
preferable. Two extreme cases, "full WW information and no fringes when measuring the 
population of quantum states" and "perfect fringe visibility and no WW information" have 
been clarified in textbooks and demonstrated with many different kinds of quantum ob- 
jects including photons [4], electrons [5], neutrons [6], atoms [7] and nuclear spins in a bulk 
ensemble with NMR techniques [8]. In Ref.[8] we further proved theoretically and experi- 
mentally that full WW information is exclusive with population fringes but compatible with 
coherence patterns. 

In order to describe the duality in the intermediate regime "partial fringe visibility and 
partial WW information" , quantitative measures for both the fringe visibility V and WW 
information are required. The definition of the former is the usual one. In variants of two-slit 
experiments different WW detectors or markers, such as microscopic slit and micromaser, 
are used to label the way along which the quantum entity evolves. A quantitative approach 
to WW knowledge was first given by Wootters and Zurek [11], and then by Bartell [12]. 
Some relevant inequalities to quantify the interferometric duality can be found in a number 
of other publications [2,13-16]. Among them, Englert [2] presented definitions of the pre- 
dictability P and the distinguishability D to quantify how much WW information is stored 
in the marker, and derived an inequality D 2 + V 2 < 1 at the intermediate stage which puts 
a bound on D when given a certain fringe visibility V . Although the quantitative aspects 
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of the interferometric complementarity have been discussed by a number of theoretical pa- 
pers, there are just a few experimental studies, i.e., the neutron experiments [17,18], the 
photon experiments [19,20] and the atom interferometer [21]. Recently, a complementarity 
experiment with an interferometer at the quantum-classical boundary [22] was also testified. 

In this paper, we experimentally investigate the interferometric complementarity of the 
ensemble-averaged spin states of one of two kinds of nuclei in NMR sample molecules for 
the intermediate situations. We follow our approach detailed in Ref. [8] but use two non- 
orthogonal spin states of another nuclei in the sample molecules as the path markers. By 
entangling the observed spin with the marker one, interference is destroyed because it is 
in principle possible to determine the states the observed spin possesses by performing a 
suitable measurement of the marker one [2] . However, in this paper, an alternative geometric 
strategy of measuring D is given and the relationship between D and the entanglement of 
the spin states is clarified. And finally the duality relation D 2 + V 2 = 1 for various values 
of D and V is testified. 



II. SCHEME AND DEFINITION 



Our experimental scheme can be illustrated by a Mach-Zehnder interferometer (shown in 
Fig. 1), a modified version of the two-slit experiment. The observed and marker quantum 
objects, represented by B and A respectively, compose a bipartite quantum system BA. 
Suppose the input state of BA to be \ipo) = |0) B \0) A = 1 00) , with |0) being one of two 
orthonormal basis |0) and |1) of B and A. Firstly, a beam splitter (BS) splits |0) B into 
+ I-Ob) ' meaning that the observed system B evolves along two paths |0) B and 
|1) B simultaneously with equal probabilities. In the meantime, path markers (PM) label 
the different paths |0) B and |1) B with the marker states \m + ) A and |m_) A correspondingly. 
The joint action of the BS and PM denoted by operation Ui, thus transforms j^o) into 

= -]=(\0)B\m + ) A + |l) B |m_) A ). (1) 
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Secondly, phase shifters (PS) add a relative phase difference between the two paths, which 
are then combined into the output state \ip 2 ) by a beam merge (BM). The joint action of 
the PS and BM, which is applied on B solely, is accomplished by a unitary operation 



1 e v 



y-e"^ f j 
And the output state \ip 2 ) = C^l^i) could be read as 



(2) 



l^2> = \ [|0) fl (|m + ) A + e^|m_) A ) + |1) B (|m_) A - e^|m + ) A )] . (3) 

Finally, measuring the population, /, of B in the state \0)b and \1)b gives 

I(<j>) = ±(l±Be( A {m + \m-) A )e i +). (4) 

where "±" correspond to the population in \0) B and \1)b, respectively. Repeating the 
measurements at different might produce population fringes. Suppose the marker states 
\m±)A = cosip±\0)a + sirup ±\1)ai and from the usual definition of the fringe visibility V = 
(I max - I min )/ (I max + I min ) and Eq. (4) one gets 

V = \ A {m + \mJ) A \ = |cos ip\ , (5) 

where <p — ip_ — <p + . 

Englert [2] proposed a quantitative measure for D by introducing a physical quantity 
Lw — the "likelihood for guessing the right way" , which depends on the choice of an observ- 
able W, 

L W = Y, max {p (W h \0) B ) , p (W h \1) B )} , (6) 

i 

where p(Wi, \0) B ) and p(Wi, \1) B ) denote the joint probabilities that the eigenvalue Wi of 
W is found and the observed object takes path |0) B or \1) B . For example, for the state 
of Eq. (f), an optimal observable W opt can be found to maximize Lw = (f + |sin</?|) /2 
in the experiments [21] and by the definition of the distinguishability D of paths D = 
— 1 + 2 m&x w {L w } [2], one gets 



D {if) = |sin</?| . (7) 

Here, we present an expression for D in an intuitively geometric way. To this end, one 
projects the marker states \m±)A into an appropriate orthonormal basis {\/3 + )a, \/3-)a}, 

\m+) A =7+|/3+)a + 7-I/5-)a, 
\m-) A = 5 + \(3 + }a + S-\(3-) a , 

where |7+| 2 + |7-| 2 = |<5+| 2 + |<5_| 2 = 1. In the two-path case the criterion of choosing 
{\P + )a, \P-)a} is to make the difference of probabilities of measuring the two states \m + )A 
and \m-) A on the basis \/3+)a to be equal to that while measuring \m + )A and \m-) A on 
\P-)a- These probability differences are then defined as the distinguishability 

L>=||7 + | 2 -|<M 2 | = ||<5-| 2 -|7-| 2 |- ( 9 ) 
The basis {\/3 + )a, \(3-)a} can be rewritten into in the computational basis: 

\/3 + ) A = cos6\O) A + sin0\l) A , 

\f3_) A = sin9\0)A-cos9\l) A , 
where 6 is the angle of the state vector \f3 + ) A with respect to the basis \0)a- In order to 
satisfy Eq. (9), from Fig. 2 and by the geometric knowledge 9 = Lp+ + Lp - — | must be held, 
which yields 

7+ = <L =cos(f -f), 
7 _ = S+ = sin(f - f ). 

Here, — ip- — tp + is the angle between the two marker state vectors in the Hilbert space. 
So from Eqs. (9) and (11), the distinguishability is equally given by Eq. (7). It can also 
be seen that the desired basis {\(3 + ) A , \/3-)a} deduced by our geometric strategy is just the 
eigenvectors of the optimal observable W opt [21]. 

These expressions for V and D are consistent with those in Ref. [21] and lead to the 
duality relation 

D 2 (<p) + V 2 (<p) = l. (12) 



Eqs. (5) and (7) reveal the sinusoidal and cosinusoidal behaviors of D and V, respectively, 
on the angle p> between \m + )A and \mJ) A in the Hilbert space H A . D and V, therefore, are 
determined by the feature of \m + ) A and \m-) A , especially by the value of ip. However, for 
any value of ip the duality relation (13) holds when two evolution paths, \0)b and \1)b, are 
labeled by quantum pure state \m + ) A and \mJ) A . Generally Eq. (12) should be replaced by 
D 2 + V 2 < 1 [2,15,16]. 

As WW information of the observed system B is stored in the states of the marker sys- 
tem A through the interaction and correlation of A and B, the distinguishability of the B's 
paths depends on the feature of the marker states, or more exactly, the correlation property 
of the combined system AB. It would be natural to examine the relationship between the 
entanglement of the system AB and the distinguishability. For a bipartite pure state the en- 
tanglement E can be denoted by the von Neumann entropy S [24], S = S (p^) = S (p( s >) , 
with S (p A ^) = -Tr(p A ^ \og 2 p A(B ^) and p A ^ = Tr A(B) (p AB ) for each subsystem. The 
entanglement E for the pure state [^i) shown in Eq. (1) is then derived as 

BM = -i^l„ g2 (l^)-i±^lo g2 (l±^). (13 ) 

It can be obtained from Eq. (13) that E = for ip = kir and E — 1 for (p — (2k + 1) n/2 
with k — 0, 1, 2, • • • , which correspond to D = and 1, respectively. A detailed quantitative 
analysis of E will be given later (see Fig. 3 below). 

III. EXPERIMENTAL PROCEDURE AND RESULTS 

The scheme stated above was implemented by liquid-state NMR spectroscopy with a 
two-spin sample of carbon-13 labeled chloroform 13 CHC1 3 (Cambridge Isotope Laboratories, 
Inc.). We made use of the hydrogen nucleus ( 1 H) as the marker spin A and the carbon nuclei 
( 13 C) as the observed spin B in the experiments. Spectra were recorded on a BrukerARX500 
spectrometer with a probe tuned at 125.77MHz for 13 C and at 500.13MHz for X H . The spin- 
spin coupling constant J between 13 C and X H is 214.95 Hz. The relaxation times were 
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measured to be 7\ = 4.8 sec and T 2 = 3.3 sec for the proton, and Ti = 17.2 sec and 

T 2 = 0.35 sec for carbon nuclei. 

At first, we prepared the quantum ensemble in an effective pure state po from the thermal 

equilibrium by line-selective pulses with appropriate frequencies and rotation angles and a 

magnetic gradient pulse [25]. p has the same properties and NMR experimental results as 

the pure state \ipo) = 1 00) . Then we transferred po to another state p\ equivalent to the state 

l^i) shown in Eq. (1) for accomplishing the BS and PM actions by applying a Hadamard 

\ 



transformation H B = ^ 



1 1 

'-V 



on spin B and two unitary transformations 



Pi = exp-iB^ +) 

(14) 

P 2 = exp(-iE7f<7jV_), 

where a* (77 = x,y,z) are Pauli matrices of the spin i, = |(1 2 ± a*) and 1 2 is the 
2x2 unit matrix. These operations were implemented by the NMR pulse sequence Ya(<p+ + 
iP-)X a (t;) Jab( ( P-- ( P+)X a (-^ )X b {t{)Yb(\) to be read from left to right, where Y A (ip++ifi-) 
denotes an </?++</?_ rotation about y axis on spin A and so forth, and Jab{}P-— f+) represents 
a time evolution of (</?_ — (P+)/hJab under the scalar coupling between spins A and B. 
Finally, the PS and BM operations were achieved by the transformation U 2 , which was 
realized by the NMR pulse sequence Xb {—0\)Yb (0 2 )Xb (— #i) with 6>i = tan _1 (— sin 0), 
and 9 2 = 2sin~ l (—cos(f)/\f2). 

In our experiments, two sets of experiments for a given value of ip = </?_ — </?+ 
were performed to measure the fringe visibility V and the distinguishability D. In the 
experiment of a quantitative measure for P, whether it is defined by the geometric 
way or the maximum likelihood estimation, the joint probabilities p (\fl±) A , |0) B ) and 
p(\P±)a, |1)b) must firstly be measured. We performed the joint measurements by a two- 
part procedure inspired by Brassard et al. [26]. Part one of the procedure is to rotate 
from the basis {|0)b|/3+)a, \0)b\P-)a, \1)b\/3+)a, \1)b\P-)a} into the computational basis 
{|00), 1 1 ) , 1 10) , |11)} (omitting the subscripts A and B), which was realized by the unitary 
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operation 

( 

Rr = 



cos a -sma 
sin a cos a 



(15) 



y sin a cos a j 

where a — | — , corresponding to the NMR pulse Y B (2a). Part two of the procedure is 

to perform a projective measurement in the computational basis which could be mimiced by a 
magnetic gradient pulse along z-axis [27]. Accordingly, , the joint probabilities p (\(3±)a, |0) b ) 
and p (\P±)a, were obtained with reconstructing the diagonal elements of the deviation 
density matrix by quantum state tomography [28]. The results are shown in Fig. 3. In our 
geometric strategy, it can be obtained from Eqs. (1) and (8) that, the information of 7+, 7_ 
or 8 + , 5_ are determined by the population probabilities, i.e., |7+| 2 = 2p(\0) b \(3 + )a), ll-f — 
2p(\0) B \P-) A ) and \5 + \ 2 = 2p(\l) B \(3 + ) A ) =, \5.\ 2 = 2 P (\1) b \P-)a). Finally, we used Eq. (9) 

1 17+1 — 1^+1 | + — |7-I I) /2 to give data points of D 
which shown in Fig. 4. On the other hand, utilizing data points of Fig. 3, we achieved the 
experimental values of the likelihood L w from Eq. (6) and obtained the D measure with 
the maximum likelihood estimation strategy, which is the same outcomes as that in our 
geometric strategy. Therefore, the intuitively geometric strategy gives the equally effective 
measure of the distinguishability D. 
For 

measuring V, we repeatedly applied the NMR pulse sequence X B (—9 1 )Y B (Q 2 )X B (—Q 1 ) 
that represents the U 2 (0) operation for various values of and detected the population of 
B in the state p 2 equivalent to the output state |^ 2 ) • A se t of appropriate values Q\ and 9 2 
were chosen to vary the values for from to 27r. Using the same reading-out pulses and 
tomography method as in the measurement of D, we reconstructed the populations of B for 
various values of <fi. The variation of the normalized populations versus showed a desirable 
interference fringe, from which the value of V was extracted. Care should be exercised in 
processing the spectra data of the different experimental runs in order to get the normalized 
populations of the deviation density matrix. 
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The objective of the present paper is to study the interferometric complementarity in the 
intermediate regime with two non-orthogonal marker states, so the experimental procedure 
mentioned above was repeated for different ip. Without loss of generality, we assumed ip + — | 
and changed the ip values from to 57r/4 by varying the ip- value with the increment of 7r/16. 
The measured values of V (ip) and D (ip) in two sets of independent experiments were plotted 
in Fig. 3, along with the theoretical curves of V (ip) , D ((f) and E (if) . The experimental 
data and theoretical curve for D 2 (if) + V 2 (if) were depicted in Fig. 4. 

From Figs. 3 and 4 some remarks can be made as follows. 

1) For ip = kn, (k = 0, 1,2, • • •) which means \A(rn + \m^) a\ = 1, two marker states are 
identical (differing with an irrelevant phase factor possible), and the state of the system AB 
is completely unentangled (E — 0). In this case no WW information of system B is stored 
in system A so that two evolution paths of B is indistinguishable (D = 0) and perfect fringe 
visibility is observed (V = 1). For if = (2k + l)n/2, i.e., ^(m+lm-)^ = 0, the marker states 
are orthogonal, and the state of the system AB is completely entangled (E = 1). This leads 
to full WW information (D = 1) and no interference fringes (V = 0). These two extremes 
are exactly the same examples that we have studied in Ref. [8] with a NMR bulk ensemble 
by population measurements. 

2) When if equals other values than kn and (2k + 1)tt/2, which corresponds to < 
\A{ m +\ m ~) a\ < 1 ? the marker states are partially orthogonal, and the state of the AB 
system is partially entangled (0 < E < 1). In this intermediate situations partial fringe 
visibility (0 < V < 1) and partial WW information (0 < D < 1) are resulted. Nevertheless, 
the interferometric duality still holds as in the extreme cases. 

3) In the whole range of ip, E varies synchronously with D. The reason is that the 
increase of E means more correlation between system B and A and more WW information 
of B stored in A, so D rises, and vice versa. On the contrary, the variation trend of E versus 
ip is opposite to that of V versus if. As the function of E versus ip has a complicated form 
there is no similar relation between E and V to the duality of D 2 + V 2 = 1. 

4) The measured values of V, D and thus the derived values of D 2 + V 2 are fairly in 
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agreement with the theoretical expectation. The discrepancies between the experimental 
and theoretical values of V, D and D 2 + V 2 in some data points, estimated to be less than 
±10%, are due to the inhomogeneity of the RF field and static magnetic field, imperfect 
calibration of RF pulses, and signal decaying during the experiments. 

IV. CONCLUSION 

In conclusion, we have experimentally tested the interferometric complementarity in a 
spin ensemble with NMR techniques. In addition to two extremes, the intermediate cases 
that the fringe visibility V reduces due to the increase of the storage of WW information 
are emphasized. The measured data of D and V in our NMR experiments are in consistent 
with the duality relation. In particular, the close link among D, V and the entanglement 
of the composite system consisting of the observed and marker states is explicitly revealed 
and explained. Though the experiment was not strictly limited in the one-photon-at-a- 
time fact, it was performed on a quantum ensemble whose dynamical evolution is still 
quantum mechanical. Therefore, our experiment provides a test of the duality relation in 
the intermediate situations. 
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Figure Captions 



Fig. 1 Schematic diagram of a two-way interferometer. The input, say |0)b|0)a, is split 
to two ways by the beam splitter (BS), then labeled by the path markers (PM), phase shifted 
by the phase shifters (PS) and finally recombined into the output by the beam merger (BM). 

Fig. 2 The state vectors in Hilbert space for defining D. {|0), |1)} represents the standard 
orthonormal basis. Two marker states {|m±)} and another orthonormal basis {|/3±)} are 
determined by the angle </?+,</?_ and 9, | + 9, respectively. The dashed line denotes the 
angle bisector between the two states {|m±)}. From the map, one can get the relation 

n ip++ip- _ 7r 

° ~ 2 4' 

Fig. 3 Normalized populations versus the angle ip, between two marker states in the 
experiments to measure D. Data points +, 0>* and H denote the joint probabilities 
p(\P+)a, |0)b), p(\P-)a,\ q )b), p(\P+)a,\ 1 )b) and p(\p_) A ,\l) B ), respectively. Theoretical 

I 1 2 

curves expressed by p(\P + )a, \0)b) = p(\P-)a, |1)b) = |cos(f - f)| /2 and p(\P-)a, \0)b) = 
p(\P + )a, |1)b) = | sm (f — f)| 2 /2 are depicted with the solid line and the dashdotted line, 
respectively. 

Fig. 4 Visibility V (denoted by O) and distinguishability D (denoted by *) as a function 
of (p. The solid lines are the theoretical expectations of V and D and the dashed line denotes 
E expressed by Eq. (11). 

Fig. 5 Experimental test of the duality relation based on the data from Fig. 3. D 2 + V 2 
is plotted as a function of (p. The solid line represents the theoretical expectation. 
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